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1 1. INTRODUCTION 

' 

The first Few Body (FB) conference was held in London in 1959, second in Brela 
(1967) and — after conferences in Birmingham (1969) and Budapest (1971) — Los Angeles 
conference in 1972 started a series of conferences held every 2-3 years. European FB 
conferences are held regularly since 1971. In 1975 European FB community was organized, 
followed by the FB research group of the American Physical Society and recently by 
Asian-Pacific FB conferences. Few Body Systems - a journal devoted to FB problems - 
was initiated by W. Plessas in 1986. FB conferences are attended by about 250 - 350 
participants from over 40 countries. FB research is more international involving over 
20 institutions from 10 countries (and more, e.g. P*ANDA involves 370 physicists from 
47 institutions in 16 countries Through its variety FB field constantly attracts 

young researchers. While initial conferences were devoted to nuclear and particle physics, 
today they include particle, nuclear, atomic, molecular, solid state physics, chemistry and 
astrophysics. 

Why did it all start? In 1953 Bethe estimated that more man-hours was devoted to 
solving nuclear force problem than to any other scientific problem in history |2j. Yet, in 
1960 "scarcely physics owed so little to so many" [3 . Half a century ago the Los Alamos 
group 0] and then Zagreb group [5] did many few nucleon studies culminating in the 
first determination of the neutron-neutron 1 So scattering length, a nn , from 2 if (n,p)nn 
and 3 i7(n,d)nn yielding —21 ± 1 fm and — 18 ± 3 fm [Hj, respectively, with no theoretical 
uncertainties included. FB studies yielded nuclear data necessary for nuclear energy, 
e.g. fusion [7,. Kinematically complete (exclusive) measurements by the Rice and 
BNL groups demonstrated quasifree scattering (QFS) and sequential decay - final 
state interactions (FSI). UCLA, Caltech, Berkeley and other groups studied 2 i?(p,pp)n, 
nucleon-nucleon (NN) Bremsstrahlung (NN7) and excited states of 4 He. 4ir measurements 
were initiated by the IKO group developing BOL detector system jTU] and by the Zagreb 
group using nuclear emulsions [TT] . 

Discoveries of heavier mesons in the early 60-ies gave birth to one-boson-exchange 
(OBE) potentials ^2] -The founding fathers of nuclear physics were aware of the impor- 
tance of three nucleon forces (3NF): "the description of nuclei in terms of NN force except 
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for the deuteron cannot be considered satisfactory" . The number of m-body interactions 
in an n-body system is n!/m!(n-m)! and that is already for 6 Li larger than n!/2!(n-2)! 
and it grows with A ^3]. I n 1957 the first 3NF was developed [T4] . 

Quantum three body problem was solved exactly by L. Faddeev in 1960 [15J followed by 
numerical results by Mitra [THj and Amado [T7] , albeit using S-wave separable potentials. 
Based on Faddeev theory Sandhas et al [THj in 1967 developed the AGS equation - our 
basis ever since. 

During this half a century a remarkable progress was achieved in nuclear forces, in 
three body calculations and in experiments. New ideas, new approaches and surprises 
marked this period extending the breadth of the FB field. I will attempt to give glimpses 
of this remarkable story which is being told so eloquently by each of you by your own 
contributions to this successful endeavor. 

2. PARADIGM CHANGE IN THE NUCLEON-NUCLEON STUDIES 

Current NN potentials form two groups, see Table 1. 1) High precision potentials |20|22[ 
123126125] fit about 6000 NN data up to 350 MeV with a x 2 / datum about 1-0, comparable 
to the phase shift analysis (PSA) [19J. These potentials need 30-40 parameters. The 
CD-Bonn potential, based on OBE model, is charge dependent (CD) and it reproduces 
pp data below 350 MeV with % 2 /datum = 1.01. AV18 includes electromagnetic (EM) 
forces, one-pion-exchange (OPE) and short term phenomenology. EM includes one gamma 
and two gamma exchanges, vacuum polarization, Darwin-Foldy and magnetic interaction 
(MI). All can be written as a sum of 18 operators (thus AV18). Potentials from inverse 
scattering (JISP) [2U| adjust their off-shell behavior to mimic 3NF. JISP6 and JISP16, 
where numbers 6 and 16 indicate whether off shell has been fixed by fitting binding 
energies (BE) of A=6 or A=16 nuclei, lead to rapid convergence in the No-core Shell 
Model (NCSM) |21| . The composite structure of hadrons and relativistic dynamics lead 
to nonlocalities at short range and Doleschall's potential is non-local at r<3 fm with OPE 
tail [22] • Kukulin's potential j2H| is a realization of a hybrid quark-meson model. Several 
3NF are used: Tucson-Melbourne (TM) [2Z|> TM99' [2E| improved to be consistent with 
chiral symmetry, Urbana/Illinois: UIX [2H] and IL1-5 [3U]. An alternative mechanism 
of generating 3NF uses explicitly A-isobar in a coupled-channel approach |3T] . The force 
employed in the isobar approach, referred as CD-Bonn + A, is purely nucleonic in T=0 
states and the coupled-channel two-baryon acts in T=l states where few constants in 
CD Bonn are retuned. It reproduces NN data with x 2 /datum of 1.02. 2) NN, 3N, 
4N.. potentials based on chiral perturbation theory (%PT). The QCD Langrangian L 
for massless u and d is chirally symmetric. The effective L starts with the most general 
L consistent with all symmetries particularly broken chiral symmetry of QCD. At low 
energies effective degrees of freedom are nucleons and pions, while heavy mesons and deltas 
are integrated out. In the chiral effective field theory (EFT) the contributions are made 
up by a series in (Q/A) n , Q= momentum or pion mass, A = 1 GeV, n>0. For a given n 
the number of contributing NN, 3N, 4NF diagrams is finite and calculable. Uncertainty at 
order n can be estimated by calculating the next order. Parameters in %PT are constrained 
by 7rN and scatterings and consistent among NN, 3N, 4N... Table 1 shows that N3LO 
is needed to get a good x 2 /datum. ;\PT generates NN, 3N, 4N etc from the same effective 
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L. The 3NF appears for the first time in N2L0 explaining why 3NF is weaker than the 
NN force. Due to partial cancellations 3NF becomes noticeable, e.g. in BE of 3 H 0.5-1.0 
MeV out of 8.5 MeV, of 6 Li 6 out of 32 MeV and in scattering processes, see V.2 and 
VII. The 3NF that appears in N2LO differs from UIX and IL 3NF, which consists of 
chiral terms of higher order. xPT is fundamental and quantitative and its appropriate 
use is equivalent to using QCD. Though %PT can be carried to an arbitrarily high order, 
predictive power is limited due to increase in the number of parameters |34|33|32] . The 
excellent description of NN data and the promising procedure (%PT EFT) are the work 
of few groups led by J.J. de Swart, R. Machleidt, V. Pandharipande, P.U. Sauer and U-G. 
Meissner. 

3. PARADIGM CHANGE IN FEW BODY THEORIES 

Only approximate calculations of the three body equation existed until Glckle and col- 
laborators achieved a breakthrough - rigorous three body calculation (r3Bc) [35] employing 
high precision NN forces, including TM, TM99', UIX and ILi, and now the Coulomb force 
|3b|37|38|39| and NnLO forces. Benchmark studies proved that various r3Bc are consistent 
within 1%. The r3Bc in most cases predicts the data very well. Of course, EM forces have 
to be included also at high energies (130 MeV) when breakup configurations include small 
angles (15°) 40 . We can use r3Bc in planning experiments. Extensive studies at 13, 65, 
135 and 200 MeV on sensitivity to NN and 3NF were done [33] . Sensitivity increases with 
incident energy. Some observables are very sensitive, e.g. at 200 MeV tensor analyzing 
powers up to 400%. These areas of high sensitivity are concentrated in specific regions 
and it is a credit to early researchers who intuitively or experimentally established them, 
e.g. symmetric constant relative energies (SCRE) (among them space star), collinearity 
|4"2"] . FSI and QFS (in early 70-ies Tjon [3H] studied sensitivities with S-wave potential). 
Several different groups are doing r3Bc |44|45|22|46j . 

FB studies include also systems with four and more bodies. Faddeev theory was ex- 
tended by Yakubovsky |4"7| . Sandhas' group |4"S] . Fonseca |39] . Pisa group [SU] and by 
developing NCSM and Green Function Monte Carlo (GFMC) [S3]- Richness and com- 
plexity increases as A increases. A=4 exhibit BE saturation, excited states and rearrange- 
ment processes. Hep process is important for solar models and Big Bang Nucleosynthesis 
(BBN)|H2]- Results of several calculations based on Faddeev theory, GFMC and NCSM 
are listed in Table 2 and a remarkable agreement is obvious. AV18+UIX still underbinds 
nuclei starting with 6 He and the discrepancy increases with A. Fitting energies of 17 
states in A<8 new 3NF were constructed: IL 1-5 which contains besides 2PE also 3PE 
and As. 3NF potentials IL1-5 are CD. Comparison with 17 states and AV18 gives aver- 
age deviation of 7.32 MeV. When the 3NF is included, the average deviation decreases 
to 2.02 MeV for AV18+UIX and to 0.04 MeV for AV18+IL3. Most significant is 10 B, 
where AV18 predicted two 1 + states below the experimental 3 + state, AV18+IL2 predicts 
correct ordering and correct energies. Similarly, for 9 Be IL2 is causing an inversion of 
levels bringing them into agreement with data. This is due to an increased spin-orbit in 
IL2. The NCSM also gives correct level ordering for 10 ' n > 12 B and 12 N when the 3NF is 
included [53. . 3B theory was applied to complex nuclear processes [S3] and to hypernuclei 
e.g. to low lying levels of \ Be treating it as the 3B system: Aaa . 
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4. TECHNOLOGY WE DEVELOP CHANGES OUR RESEARCH AND US 

Developing and building instruments played an important role throughout the history 
of science. During this half a century major developments occurred in: accelerators, 
detectors, electronics, computers and data analysis. Few highlights follow. Review of 
facilities in Europe, the USA, Russia and Japan is given in refj^S]- There are now 17,500 
accelerators: 120 producing beams of energies larger than 1 GeV, about 1,000 low energy 
accelerators for research, 100 synchrotron radiation facilities, more than 7,500 radiother- 
apy accelerators (mainly linacs) and more than 7,000 for ion implantation and condensed 
matter work jH7j. Tandems played important role in our research in 1960-90- ies and few, 
notably TUNL and Cologne, are still leading centers. Early SF cyclotrons, then LAMPF, 
TRIUMF, PSI, IUCF, RCNP, MIT/BATES and now JLab (to be upgraded to 12 GeV), 
KVI, HI7S (High Intensity Gamma Ray Source at Duke), MAMI -Mainz Microtron to 
be upgraded to 1.5 GeV, CELSIUS - accelerator and storage ring at Uppsala, LEGS - 
Laser electron gamma source of highly polarized monochromatic gamma 150 - 470 MeV 
by Compton backscattering from 2.8 GeV electron beam at National Synchrotron Light 
Source at BNL, ESRF - European Synchrotron Radiation Facility, storage ring at Upp- 
sala and electron storage ring TERAS at AIST, Tokyo are playing a crucial role in our 
research. Radioactive ions beams (RIB) opened new insights. For the first time it is 
possible to use beams of N/Z = 3, e.g. 8 He facilitating the transfer of 4n. RIB are pro- 
duced by: i) isotope separation on line using two independent accelerators: accelerator 
for production of short lived radio-nuclei in a thick target connected to an ion source 
and a post accelerator, ii) in flight which relies on very energetic beams of heavy ions 
impinging on a thin target and resulting in fission or fragmentation with velocities of 
outgoing particles comparable to velocity of the incident projectile. Cocktail of different 
particles identified by mass, charge and momentum are produced and do not need further 
acceleration. And iii) a combination of i) and ii). A variety of RIB are available with very 
short lifetimes, e.g. 14 Be (T= 4.2ms), 8 He (T= 119ms), 6 He (T= 0.807s), 8 Li (T=0.838s) 
and long-lived 10 C (T= 19.3s) and 18 Ne (T=17s). High beam intensities are achieved: up 
to billion particles per second. The number of RIB facilities is increasing. Not a complete 
list includes ISOLDE at CERN, SPIRAL at Ganil, FAIR (construction begins in 2008) at 
GSI, ISAC at TRIUMF, RIBF at RIKEN (under construction) and SBL (both in Japan), 
NSCL at MSU, Louvain-la-neuve, Orsay, ORNL, Berkeley, FSU, RIBRAS (Sao Paolo), 
RIBLL (Lanzhou, China), Argonne and Catania [58J. DOE is considering preliminary 
engineering design of advanced exotic beam facility in 2011. Several kaon beams were 
developed strengthening hypernuclei research: e.g. COSY, JLab, Nuclotron at Dubna 
and DA$NE (Double Annual ring For Nice Experiments), Frascati, colliding e~ and e + 
beams to produce $(1020) decaying into neutral and charged kaons. DA$NE is delivering 
about 150 $/s. Beams of antiparticles, e.g. antiprotons and antideuterons are developed. 

Detectors range from single to complex structures matching in R&D and cost those of 
accelerators. Detectors are used in stand-alone mode as in underground neutrino experi- 
ments or with accelerators. A stand-alone detector is the proposed Majorana 76 Ge double- 
beta decay detector: a 500 kg of 76 Ge isotopically enriched to 86% allowing a sensitivity of 
half-life measurement of 76 Ge of 4xl0 27 y over a data acquisition period of 5000 kg years 
[55] . Some examples of detectors used with accelerators are: JLab CLAS, SALAD-KVI 
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(small angle large acceptance detector) (AGOR facility), now BINA, RIKEN spectrograph 
SMART; Gammasphere array of 100 large Compton-suppressed Ge crystals, GRETA and 
GRETINA (gamma ray energy tracking in-beam nuclear array), Crystal Ball (at AGS and 
now at MAMI), Crystal Barrel and TAPS at ELSA, WASA (wide angle shower appara- 
tus) at CELSIUS and multidetector systems developed at Erlangen, TUNL and Cologne, 
superconducting magnetic spectrograph S800 at MSU, Blowfish detector array consisting 
of 88 neutron detectors, ANKE detection system consisting of range telescopes, scintilla- 
tion counters and MWPC at COSY and detection systems developed at JLab. Detection 
systems for hypernuclear research are hyperball, FINUDA and P*ANDA. 

Astronomical observations require measurements of high energy gammas up to 10 20 eV, 
rapidly varying gammas in TeV range, and neutrinos. Underground neutrino detection 
system include Kamiokande, SuperKamiokande, BAKSAN, SNO, KamLAND and a next 
generation of solar neutrino detectors BOREXINO in Gran Sasso, detectors under the 
South Pole (Amanda and Ice Cube), underwater detector projects in the Mediterranean 
(off the coasts of France and Greece) and in Russia. 

Computers - for high speed complex computing, for automated measurement and con- 
trol and for storage, retrieval and data analyses - coupled with networking, internet, www 
and fusing computers and telecommunications - revolutionized research. Higher energies 
and higher intensities of accelerators open the possibility for processes generating vast 
amount of information that has to be stored and analyzed and total accumulated data in 
a typical experiment are now orders of magnitude larger than few decades ago. Software 
cuts on kinematical conditions and particles ID reduce the data stored on tape by several 
orders of magnitude. Many studies would be impossible without Monte Carlo simulations. 
The progress from tube computers that were used in the 50-ies with punched cards to 
present supercomputer, DNA computer constructed in 2004 by E. Shapiro at Weizmann 
Institute and a prospect for quantum computing made a paradigmatic change. 

5. FEW BODY RESEARCH - RESULTS AND PROBLEMS 

There is an impressive amount of FB data. NN and nucleon-deuteron (Nd): proton- 
deuteron (pd) and neutron-deuteron (nd) elastic and breakup data: total, partial and 
differential cross sections, spin observables measurements involving N and d. The long- 
standing issue of a nn will be addressed by doing nn scattering experiment (first tried with 
underground nuclear explosion) using the pulsed aperiodic Yaguar reactor at Snezhinsk 
giving neutron flux of density of over 10 18 cm -2 s _1 jHO]- Some results and problems: 

V.l. Symmetries: 

V.l.l. Unitarity of the Cabibbo - Kobayashi - Maskawa matrix The Particle Data 
Group reported in 2002 V u d = 0.9734 ± 0.0008 (from beta and neutron decays), V us = 
0.2196 ±0.0023 (from K+ and K° decay) and V ub = 0.0036 ± 0.0007 (from semileptonic B 
decay), wherefrom breaking of unitarity occurred: Delta = 1 — (|U u d| 2 + \V US \ 2 + \V u b\ 2 ) = 
0.0043 ±0.0019, i.e. 2.2 a |B2|. Several measurements done after 2002: BNL E865, NA 48, 
KLOE, KTev and ISTRA± including K and hyperon decays gave V us = 0.2261 ± 0.0021 
leading to Delta = 0.0004 ± 0.0011 compatible with unitarity [03] • Studies of pion decay 
and of 12 superallowed transitions + — ► + from 10 C to 8 Rb are currently performed 
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V.1.2. CPT invariance jHIj requires the equality of masses and lifetimes of particles 
and antiparticles. The upper limit of 10~ 18 is established based on the equality of masses 
of K° and its antiparticle. 

V.1.3. CP violation [HI] in the quark sector is known since 1964 and research goes on 
(evidence in — > Knn; permanent electric dipole moment (EDM) of any fundamental 
particle violates P and T, limits on EDM for electron is 1.6 x 10 -27 ecm, for neutron 
6.3 x 10~ 26 ecm and for proton — 3.7±6.3 x 10 -23 ). However, nothing is known about 
CP violation in the lepton sector. Based on solar and atmospheric neutrino and reactor 
antineutrino measurements neutrinos do have mass (otherwise there would be no oscilla- 
tion, albeit according to WMAP observations the sum of masses of three neutrino flavors 
is less than 1 eV [HI])- Comparison of behavior of neutrinos and antineutrino provides 
information on CP violation in the lepton sector. 

V.1.4. Conservation of lepton number: neutrinoless double beta decay (Z,A) — > (Z + 
2, A) + e~ + e~ studies yield for the upper value of the 76 Ge lifetime 1.9 xlO 25 y with a 
confidence limit of 90% [61 J. 

V.1.5. Time variation of the fine structure constant olqed'- observation of the variability 
of basic constants could give information on the existence of extra dimensions jHS] • BBN 
offers opportunity to test time variation of fundamental constants since reaction rates 
depend on fundamental constants [UHl- Crucial rate is that for np —* at 30 - 130 keV. 
The 5% uncertainty associated with this process contributes a significant uncertainties 
in the abundances of light elements produced in the BBN. This process was measured 
at HPyS and the data jHZj are in agreement with the EFT calculation [HE] and with the 
NN potential model of ref jJjU] thereby reducing the uncertainty of 5%. Discovery that 
the quasar absorption lines in interstellar media suggest a small evolution of the olqed'- 
a/a = (6.40±1.35)xl0 _16 /year [70] is related to the predictions of string theory and 
Kaluza-Klein theory [71]. This observation is not confirmed. Laboratory experiments 
using narrow atomic hydrogen Is — > 2s transition compared to atomic clock did not see 
any change [72]. Primordial helium abundance depending on 6 fundamental constants 
and WMAP observations put a bound on Aa/a of few times 10~ 3 . Stringent bounds can 
be obtained from deuterium and lithium primordial abundances j73j . 

V.1.6. Supersymmetric particles search [HT] resulted in lower limits for their masses, 
e.g. for neutralino 46 GeV, for chargino 94 GeV, for selectron 73 GeV, for squark 259 
GeV and for gluino 195 GeV with CL of 95%. and 

V.1.7. Charge dependence and charge symmetry breaking (CSB) These were the first 
broken symmetries studied resulting in SU(2) and SU(3) and non-Abelian gauge theories. 
On the fundamental level CD and CSB are due to differences between masses and charges 
of u and d quarks [73]. The consequences on the hadronic level are: i) the difference among 
hadron masses: m^ - m u manifests itself in 7r ± - 7r°, and n - p mass differences (m n > m p 
has significant consequences for the structure of the Universe), ii) meson mixing, p-u is 
important, and iii) irreducible meson-photon exchanges. Pion mass difference accounts for 
60% of CD in 1 So jZH- Nucleon mass splitting has several consequences. The most trivial 
is the change in the kinetic energy whereby the a nn becomes more negative by 0.25 fm. 
OBE diagrams are negligibly affected, but TPE are. The p-u mixing effect was subject 
to criticism and counterarguments and this issue is still open [73]. CSB is observed in: 
i) Sa exp = cl pp — a nn = (—17.3 ± 0.3) — (—18.9 ± 0.4) = 1.6 ± 0.5 fm (Only the a nn value 
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from the 7r~d measurement is taken -18.93± 0.27 ex ±0.3^ jTHJ- Reassessing theoretical 
uncertainties |77j leads to even smaller value of ±0.2 th . Values from nd breakup are 
discussed in V.3). Various models predict for values in agreement with this 

experimental result: chiral quark model [TS] predicts 1.37 - 1.67 fm and it correctly 
predicts a np = —23.749 fm. CSB due to p-u mixing gives 1.508 fm, the same result as 
the TPE taking into account n-p mass difference. Therefore, a pp — a nn cannot distinguish 
between two hadronic mechanisms generating CSB jj^]. ii) 3 H and 3 He BE difference is 
764 keV, mainly due to EM effects and after all corrections the 65±10 keV is attributed 
to CSB [SU]. Isospin violating 3NF within the EFT contributes 5 keV to the CSB in the 
3 H- 3 He [HI]. CSB in 3 H- 3 He is due mainly to S-wave forces: CSB potential derived from 
AV18 fixed by the best fit to pp data, and then adjusting its central force in the S=0, T+l 
state to fit a nn gives for 3 H- 3 He 62.1 keV for 1 Sq and 65.1 keV when all T = 1 states are 
included (7^1 • Results for CSB due to p-u mixing and due to n-p mass difference are 60.9 
and 57.6, respectively for 1 So and 65.7 and 60 keV when all T = 1 states are included. 
For the S=l and T=l states ( 3 Pj) it is assumed that the CSB splitting is the same as 
in S=0, T=l. iii) Nolen - Schiffer anomaly - energy difference between mirror nuclei that 
cannot be accounted for by EM and nuclear models - is of the order of 200 - 400 keV and 
50% comes from L>0 NN forces [70]. Nuclear matter is sensitive to L>0 forces. CSB 
AV18 and CSB potentials due to n-p mass difference and p-u mixing predict different 3 Pj 
phases: nn and pp 3 P and 3 Pi differ by 0.5° to 1.0° and 3 P 2 by less than 0.2° 
CD in 3 Pj waves was studied through PSA [HI] and found that np and pp 3 Po at 25 Mev 
differ by 2.2°. iv) difference between n and p analyzing powers studied at 183 MeV [SB] . 
347 MeV [HH] and 477 MeV [H7J- Theory predicts data very well [88^ and p-u mixing 
gives different contributions at different energies, v) superratio of n + ix~ on 3 H and 3 He. 
vi) total cross section for D(d, a)n° measured [89J at 228.5 and 231.8 MeV at IUCF and 
found to be 12.7±2.2 pb and 15.1±3.1 pb smaller than simplified model [SU] based on 
the EFT predictions of 23.0 pb and 30.8 pb. vii) asymmetry in n + p — > d + tt° a(9) 
measured at TRIUMF: (17.2±8st±5.5sy)10 -4 [HI]. In the framework of the EFT main 
contributions are mixing and rescattering which through partial cancellation produce an 
upper value of 69xl0~4 [§2]. viii) recent theoretical analysis of the ttN data found a 
small CSB and ix) A separation energies in and jjie, which after removing Coulomb 
effects amounts to 390 keV, i.e. about five times larger than the CSB in A=3. Extraction 
of strangeness form factors, the NuTeV anomaly and hadronic corrections to g-2 are 
related to CSB 93 . Parity-violating L-R asymmetry in e-p scattering is sensitive to s 
quark contributions and measurement at JLab|94j has to be corrected by CSB. Ratio 
of charged- current and neutral-current cross sections for neutrinos and antineutrinos on 
isoscalar target is equal to 1/2 — sin 2 9w- Measurements [95J give sin 2 ^ three standard 
deviations larger that Standard Model (SM) prediction. CSB corrections account for 
about a half of the discrepancy [96J. ^7 — 2 is known to 1 in 10 6 , but it is 1 - 2.4 standard 
deviations off SM predictions and accurate CSB effects of p-u mixing account for 20% 

EJJ- 

V.2. Evidence for the 3NF: Composite structure of nucleons and xPT EFT lead to 
3NF. We summarize experimental evidence for 3NF: 

1) High precision NN forces underbind 3 H BE and the disagreement gets worse as A 
increases. There is a correlation among 3N and 4N BE |9S] recently extended to light nuclei 
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and nuclear matter [99J. Phenomenological 3NF with parameters adjusted to predict 3H 
BE can account for BE and energy levels up to A = 16 [28 30 31 j. Using the same NN 
and 3NF a fit to 4 He is obtained HUOj indicating that 4NF gives a small contribution. 
Ordering of energy levels of light nuclei, particularly of 10 B, requires 3NF and IL2 with 
increased LS gives the best description. 

2) PSA of pd elastic scattering E p = 4 - 22.70 MeV data requires 3NF to achieve 
compatibility between 3N and NN phase shift |102j . Just as AV18 cannot predict Et, 
it fails with 2 Si/ 2 , and agreement is achieved if UIX is added. Ay and iTn depends on 
4 Pj and these are related to the NN 3 Pj. The correct description of Ay depends on 
2 f3/2 — 4 P3/2 mixing €3/2-. Small changes in eyz- are required to describe iTn and Ay. 
Essential problem is that 4 Pj required by 3N PSA are incompatible with 3 Pj demanded 
by the NN PSA. 3NF affecting 4 Pi/2 and e 3 / 2 _ are required and it seems that the LS type 
3NF from refs |103ll04j fulfills this role. 

3) R3Bc with the CD Bonn (TUH] fit total nd cross section data [HJEl up to 150 MeV. 
Above this energy and up to 200 MeV inclusion of the 3NF fits the data. Relativistic 
kinematics increases the total cross section, but boost effect outweighs it j!07j . Including 
3NF and relativity still does not account for the nd total cross section at high energies. 
Data up to 200 MeV provide evidence for the 3NF. 

4) Disagreement between the data at the minimum in the elastic Nd scattering differ- 
ential cross section and the r3Bc using NN forces is decreased when the 3NF are included 
110811 09 . Measurement of the nd a{9) at 95 MeV |110j covers the full angular distribution 
(20°-160°) measuring neutrons and deuterons at forward angles using detector SCANDAL 
at Uppsala University. Data from 100°-160° show the need for the 3NF. Relativistic ef- 
fects modify a{6) in the domain around the minimum and in general increase the cross 
section therefore having a similar effect as the 3NF. 

5) Extensive study of H(d,pp)n at 130 MeV performed at KVI with SALAD system 
|lllj in 72 kinematic configurations comparing the cross section data with NN, NN + 3N, 
isobar and N2LO (with NN and 3N) predictions, ignoring EM interactions |112j . shows 
agreement between the theory and of the data, but also some disagreements. Fits to the 
data allowing for 10% uncertainty in the overall normalization is improved if 3NF are 
added. Some discrepancies still remain. 

6) High precision NN potentials could not explain polarization transfer coefficients 
measured at 22.7 MeV |112j . but a good agreement is obtained when 3NF TM' and UIX 
are included. NLO does not fit the data, but N2LO which includes NN and 3NF does. 

7) Proton-deuteron elastic scattering observables Ay, iTn and Cyy at 197 MeV incident 
proton energy cannot be explained by CD Bonn or AV18 NN potentials but fit is obtained 
adding TM 3NF [TLlj . 

V.3. Experimental evidence not compatible with current theories 

1) There are discrepancies between 7rNN coupling constants extracted from different 
observables |114j . 

2) Cross sections 3 He(7,p)d and 3 He(7,pp)n at 10.2 and 16 MeV j!15j are compared 
with AV18 and AV18+UIX. Theory predicts 10-14% larger cross section at 16 MeV and at 
10.2 MeV predicts 3 He(7,pp)n a factor of 3 too high. Neutron energy spectra measured at 
HI7S using 12.8 MeV linearly polarized gammas peak at forward angles at lower energies 
than r3Bc with CD-Bonn and Coulomb |116j . 
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3) nd and pd Ay data revealed 25 years ago difference between them and the 3B calcu- 
lation |117|118| . As the accuracy of the data improved and additional data accumulated 
and as r3Bc now including EM are made the disagreement became more convincing. In- 
clusion of MI raises the calculated Ay for pd by 4% bringing it closer to the data, but 
decreasing by 3% for nd increasing the disagreement. MI is appreciable at energies below 
5 MeV. EM plays a significant role |37|38| . but it does not solve the Ay puzzle at low 
energy: difference between the nd and pd Ay peak values and r3Bc from very low ener- 
gies to about 20 MeV |119j . Explaining these disagreements would require a long range 
3NF. Success of the LS 3NF |103j and Kukulin's potential |12Uj require further studies. 
A remedy would be CD and CSB in the 3 Pj waves, since Ay and iTn are a magnifying 
glass for the 3 Pj force. Such attempts were made unsuccessfully |121|122j . 

4) Though Ay data agree with r3Bc at incident nucleon energies of 30 MeV, the Ay 
puzzle reappears at higher energies, i.e. at 150 and 190 MeV, where data are not repro- 
duced by high precision NN + TM |123j . Disagreement is also found in deuteron tensor 
analyzing powers |4*T] . 

5) SCRE, particularly space stars investigated at incident proton energy of 10.5 - 130 
MeV and for 10.3 - 25 MeV neutron energy by several different groups. |124ll 25 126 . 
Data are consistent: at 10.3 and 13 MeV nd has 50% larger cross section than pd. At 
10.5 and at 13 MeV r3Bc even including the isobar and EM do not explain either pd or 
nd data. 3NF and MI have negligible (< 2%) effects and the cross section is dominated 
by well known S-wave. R3Bc fit well at higher energies. 

6) While r3Bc fits well np QFS in a wide range of incident energies, it does not explain 
nn and highly accurate pp QFS data. Including the Coulomb force brings r3Bc in better 
agreement with the pp QFS data, MI effects are less than < 2% and effects of the A or 
other 3NF are small. 

7) Four-body calculations with high precision NN and UIX 3NF do not fit the n - 3 H 
total cross section data [52"] . 

8) Studies of nn FSI in nd breakup using r3Bc were done at incident energies 13 - 25 
MeV. The a nn influences the absolute cross section and the shape of the FSI enhancement 
and a nn was extracted by both approaches. Kinematically incomplete study at 17.4 MeV 
|127j gave a nn = (— 16.5±0.69st± 052sy) fm. The study at 13 MeV measuring neutrons' 
time-of-flight at 9 n \=9 n 2= 20.5°, 28°, 35° and 43°, and the deposited energy of the charged 
particle gave a nn = (— 18.72±0.13stat±0.65sys) fm |128j using absolute cross section, and 
(— 18.84±0.47) fm using the shape. Earlier studies at 13 MeV |129j did not use the r3Bc, 
but if its shape is reanalyzed it yields — 17.9±0.5 fm. Studies at 16.6 MeV and 25.3 MeV 
measuring the momenta of the outgoing p at 41.2° and n at 55.5° at opposite sides of the 
beam and consequently using a thin deuterium target gave [130 a nn = (— 16.2±0.3) fm 
and (— 16.3±04) fm, respectively. This disagreement between a nn extracted at 13, 16.6, 
17.4 and 25.3 MeV is not understood. With additional neutron detectors at 55.5°, 69° 
and 83.5° at the opposite side of the beam np FSI was studied giving a np = — 23.5±0.8 
fm |131j . In a follow-up study at 25.2 MeV measuring neutron and proton at the same 
angle of 32o gave a np = -24.3±1.1 fm |130j . These comparisons validate both studies to an 
extent that r3Bc explains np FSI. However, it is known that some np FSI configurations 
are not well described by the r3Bc. Joint effort by both groups is in progress at 19 MeV 
using simultaneously both geometries: 1) n-n coincidences measured at n i=6 n2 = 35° and 
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51.7° and 2) n - p coincidences at 9 P = 45° and 9 n = 51.7°. Preliminary result for the first 
geometry is (-17.6±0.2st±0.9syst) fm |132| . Discrepancy between results for a nn could be 
due to the inadequacy of the included 3NF. This is an old idea introduced 25 years ago 
|133| and a simple model with the Fujita-Miyazawa 3NF explained the difference between 
a nn obtained from two different configurations. Later study |134j reanalyzed data used 
in extracting a nn and found a much wider spread in a nn values and r3Bc showed that at 
specific production angle at each incident energy the 3NF effect vanish (for 13 MeV it 
vanishes around 45°). The data |131j confirm small sensitivity to the 3NF. The difference 
in a nn cannot be due to MI as suggested |135j . since r3Bc showed that they are smaller 
than 0.2 fm |136j . Namely, MI influence both NN and 3N systems. 

9) Cross section and tensor analyzing powers Ayy at 19 MeV for the reaction H(d,pp)n 
in SCRE |137j - except for the coplanar star - do not agree with the predictions of r3Bc 
using TM, TM' and UIX or the coupled channel CD-Bonn + A and including Coulomb 
forces. Calculations were also done for N2LO including 3NF. 

10) Radiative pd capture cross section, Axx and Ayy using polarized deuteron beams 
from RNCP at 140 and 200 MeV could not be explained dHSj by r3Bc. 

11) High precision NN and 3NF TM' and UIX could not explain 108 - 190 MeV a{9) 
data |139j . 

12) Disagreement between a (9), Ay and polarization transfer coefficients Kyy' and Kxx' 
in the d-p elastic scattering at 250 MeV and all high precision NN + TM. The inclusion 
of TM 3NF deteriorates the fit to polarization transfer coefficients Kxz' and Kzx'. This 
is a first complete set of pd at intermediate energies |140j . 

13) While most of the data from the reaction H(d,p)d at 135 MeV/nucleon |109j are 
explained well by including the 3NF, Axx, Ayy, Axz analyzing powers and Kyy' and 
Kxzy' could not be indicating inadequacy of the current 3NF. 

14) Two-body 3 He and 4 He photodisintegration data measured at JLab at 0.35 - 1.5 
GeV do not agree with the best available calculations, albeit these calculations are rather 
old-fashioned j!4f] . 

15) Lithium problem - a statistically significant discrepancy between standard BBN 
predictions for 7 Li and twice smaller observational value over a wide range of metallicities 
|142|143j . Existence of the metastable (T > 1000 s) negatively charged electroweak scale 
particle X~ would alter the prediction of 7 Li and other primordial abundances for A>5 
via the formation of the bound state with nuclei during the BBN. We might be tempted 
to repeat the famous statement that these discrepancies are 15 minor clouds on a bright 
sky, but we should know better. 

6. IS THERE AN END IN SIGHT FOR FEW BODY RESEARCH? 

Successes described in II and III point toward a Standard model of FB theory. However, 
FB field is broader than few nucleon and there are many open issues: 

1) A hypernuclear event caused by cosmic rays was observed in 1953 |144j . In 70-ies 
kaon beams were developed. There are 35 known A hypernuclei from A H to A 09 Bi mostly 
below A = 20. Only three double-lamb da- hypernuclei AA He, aa-^g and AA B are known 
[145J. Evidence for AA H was suggested |146j . AA hypernuclei provide information on AA 
interaction, properties of multi-strange systems and neutron stars |147j and give insight 
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into the onset AS hypernuclei stability [148J. The AA BE Baa was determined from 
AA He suggesting that effective AA interaction is considerably weaker than thought earlier. 
Baryon-baryon interactions for N, A, £ and S were studied by the Nijmegen group |149j 
and in chiral quark model |150j . Double A hypernuclei were studied in the relativistic 
mean-field theory and in the 4B model: x +a +A +A, where x is n, p, d, t, 3 He or a. 

2) Renaissance of QCD spectroscopy. Proton energy distribution from the reaction 
4He(K-stopped,p) shows a peak interpreted as a neutral 3 baryon state, M = 3117, T < 
21 MeV with isospin 1 and strangeness -1 and main decay mode NN£ [156j . The nature 
of A (1405) is not fully understood yet. It could be a 3q, a molecular-like meson-baryon 
bound state or exotic |151j . The £(1480) is described as a bump by PDG |152j . The 
Crystal Ball study of the reaction K~p — > 7r°7r°A do not show tt°A resonance, but are 
dominated by £°(1385) |153j . Recent study of the reaction pp — > pKirX, where X is 
the unidentified residue, is interpreted as evidence for a neutral hyperon resonance Y°* 
decaying into tt+X" and 7r"X + with M = 1480 ±15 MeV, T = 60±15 MeV and since it 
is neutral it can be either £° or A. Light baryonic states were reviewed and search for 
A*, £*, S* and Q* states - expected but yet undiscovered - was stressed j!54j . Relatively 
narrow deeply bound antikaon nuclear states in FB systems were suggested |155j and 
searched for |156j . The dominant problem is missing resonances - the link between genuine 
3q resonances and scattering data |157j . Agreement on confidence level of important 
resonant states is still not reached. It is suspected that missing states couple weekly 
to the formation channels (I58ll59j and it seems that it is related to the importance of 
inelastic channels as recently confirmed by the Zagreb group demonstrating that inelastic 
channels, particularly ttN — > r]N are decisive in unambiguously confirming |160| one of 
controversial states: AT(1710)Pn. Inelastic channel are very important and experimental 
efforts should be strengthened (proposals EPECUR and MIPP [161J). Mesons, baryons, 
resonances etc are topics of several conference series MENU, ETAMESON, NSTAR [162J 
and of PSA workshops recently held in Abilene, Zagreb and Tuzla. 

3) r\ physics includes studies of 77-mesic nuclei and of rare decays [163 . Charge con- 
jugation invariance (CCI) forbids the decay of pseudoscalar mesons into an odd number 
of 7s and search for 77 — > 37 is a direct test of CCI for EM interaction of hadrons. The 
upper limit for the branching ratio (BR) for this process was established by two groups 
jTMj yielding BR < 4xl0~ 5 (Crystal Ball) and BR < 1.6xl0" 5 (KLOE at DA$NE). The 
process 77 — > 7r°77 is a test for x?T. Since there is no direct order 771" coupling and the 
second order involves G-parity violating transition, the first contribution comes from the 
third order. In 1982 GAMS-2000 gave BR = (9.5±2.3) 10" 3 , and if confirmed would be 
the first failure of x?T. Recent measurement gave BR = (3.2 ± 0.9) 10 -4 , converted to T 
= (0.45 ± 0.09 si ± 0.08 sy ) eV vs. the xPT result of 0.42± 0.20 eV. CBELSA TAPS 
collaboration is redoing this study with the aim of reaching a better accuracy [166J. 

4) QCD allows systems with explicit gluonic degrees of freedom: hybrids (qq*g) and 
glueballs (gg). Hybrid configurations also appear as intermediate states |167| . Smoking 
gun for hybrids is to have JPC forbidden by qq* models, e.g. , 1~ + , 2 + ~ [168 . 
Evidence for several hybrid candidates are found: i) 1~ + in the reaction pn~ — > 7t~7t~tt + p 
at 18 GeV/c in addition to expected states 1++ al(1260), 2++ a2(1320) and 2~+ tt 2 (1620) 
Pig , ii) states at M = 1709 and 2001 MeV and at 2096 MeV from the reaction pi\- -> 
7r~7T~7r + 77p |170j and iii) exotic meson decays to uj n° n~ |171j . Evidence for glueballs are 
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not found [ 172] . Strangeness in the proton and in N(1535) suggesting an important ss* 
component, relevant e.g. in the search for Cold Matter |173| . is discussed in |174|175| . 

5) There are speculations that the short range repulsion of the NN force is increased 
and the intermediate range attraction decreased in the medium |176|177j . 

6) Halo structure |178ll79|18Uj (e.g. 14 Be |181j ). Borromean, tango and samba nuclei 
and quantum proximity resonances |185j are related to Thomas |182j and Efimov |183j 
effects. Borromean nucleus |184j is a 3B bound state where none of its 2B subsystems 
is bound, e.g. 6 He. Excited state of 12 C is a Borromean excited state of paramount 
importance in SINS. Tango nuclei have only one bound subsystem and in samba nuclei 
two of the 2B subsystems are bound. 10 Be can be treated as a 3B ( 8 Be + n +n) or as a 
4B problems, and there could be a molecular <x2n:a band |186j . Using QFS and inverse 
kinematics to study the reaction H( 6 He, ab) and H( 8 He, ab), where a and b stand for 
pn, pa and p- 6 He, momentum distributions inside 6 He and 8 He were measured |187j . 3B 
calculations are consistent with data on nn rms distances in 6 He, n Li and 14 Be [188J. 
QF reactions (A(a+b) + B -»• C +D +b (initiated in 1967 by M. Furic [MI]) and two 
spectators QFS: (A(a+b) + B(c+d) — > C + D + bs + ds) (initiated by D. Miljanic in 
1974), later studied at NRL |19U|191j . could be relevant for productions of neutron- and 
proton-rich nuclei. Studies of xH and xn are initiated with 14 Be and 8 He beams. Resonant 
states in 7 H are found by H(8He,pp) 7 H, evidence that 8 He contains 6 He subsystem in an 
excited state 2 + with a large fraction [192 and 7 He levels were studied [193J. Similar 
studies are done for hypernuclei, e.g 7 AA H |194j . 

7) Session were devoted to atomic, molecular and condensed matter physics and chem- 
istry at FB meetings starting from Nanning in 1985 |195j . FB18 included talks on pairing 
of atomic fermions [196J, FB effects in cold atoms |197j and evidence for Efimov states in 
ultracold gas of Cs atoms [198 . Chemical potentials of the two spin states are equal in 
superconductors. Mismatched chemical potentials are important in magnetized supercon- 
ductors and cold dense quark matter at the core of neutron stars. Pairing in a polarized 
gas of 6 Li atoms was observed |199j . Observation of Efimov states in nuclear systems is 
hampered by Coulomb forces and for A>4 there are no Efimov states 20 U|2Ulj . In molec- 
ular physics helium trimer is predicted to have an excited Efimov state but its existence 
is not confirmed |2U2|2U3j . Studies of 3B recombination processes in ultracold gases with 
magnetically tunable 2B interaction of Cs atoms (thus varying 2B scattering length) es- 
tablished evidence for Efimov states |2U4j . The e~-e~ correlation in initial and final states 
of double ionization and (e,3e) processes is an important tool in atomic physics [205] and 
it persists even at high energy impacts as demonstrated by the asymptotic formula for 
the ratio of double to single cross section in Compton scattering 206J. 

7. CHALLENGES 

Surveys presented in V and VI show the need for extensive experimental studies - high 
accuracy and innovative. Physics depends on measurements and several proposals for new 
major facilities and upgrading of existing facilities were submitted. It is in decisions to 
build such facilities that planning of science comes as shown in a famous Bernal - Polanyi 
debate in the 40-ties. This is contained in a prescription for success: get the best people, 
give them best working conditions and do not disturb them. 
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Table 1 



X 2 /datum for various NN potentials. 



YEAR 


N° of data Nijmegen PSA 


CD-Bonn 


AV18 


N3LO 


N2LO 


NLO 


1992 


pp(1787) 1.00 


1.00 


1.10 








1999 


pp(2932) 1.09 


1.01 


1.35 








1999 


NN(5990) 1.04 


1.02 


1.21 








0-290 MeV NN data 


1999 


np 




1.04 


1.10 


10.1 


36.2 


1999 


nn 




1.38 


1.50 


35.4 


80.1 


N° of parameters 35 


38 




24 


9 





Table 2 

Binding energies (MeV) of 3 < A < 16 calcualted with 3NF(37,37a) 



3 H 


3 He 


4 He 


6 He 


6 Li 


7 Li 


8 He 


3 H 


8 Be 


16Q 


Experiment 


8.488 


7.72 


28.296 


29.27 


31.995 


39.2 


31.598 


56.5 


127.619 


AVI 8 
F:CDB+TM 
AV18+UIX 


7.628 
8.478 
8.478 


6.917 
7.735 
7.733 


24.07 
29.25 
28.34 




26.09 










G:AV18+UIX 
AV18+IL2 


8.487 
8.43 


7.739 
7.67 


28.33 
28.27 


28.1 
29.4 


31.1 
32.3 




27.2 
31.3 


54.4 
56.6 




N:AV8'+TM' 
JISP16 


8.496 


7.797 


28.374 


28.189 
28.32 


31.00 






53.3 


133.8 


N:N2LO NN only 
Full 












34.6 
36.7-38 









Though most of discrepancies listed in V.3 were confirmed by more than one measure- 
ment, it is necessary to perform a systematic consistency check and pruning of FB data 
as was done for NN data paying attention to normalization errors. (There is currently 
a disagreement between KVI and RIKEN high energy Nd data [207 208] and both data 
disagree with r3Bc.) Some data show clear evidence for 3NF and some agree better with 
no 3NF and with no Coulomb, though are affected by both. Likely, disagreements listed 
in V.3. do not imply any new physics, but demand careful inclusion of relativity, EM, 
3NF, CD and CSB, possibly fine tuning of the NN force. Ay puzzle, n- 3 H and ordering 
of levels of light nuclei suggest the need to fine tune P waves. High precision potentials 
from Table 1 are not phase equivalent, they predict different E t and have different P^ : 
CDBonn predicts 8.013 MeV and 4.85%, respectively, AV18 7.628 MeV and 5.76%, Do- 
leschall 8.48 MeV and 3.6% and N3LO 7.855 MeV and 4.51%. They have different short 
range interactions and differ conceptually. Moscow- Tbingen approach [209J introduces a 
dibaryon and claims that dibaryon degrees of freedom produce effects in FB systems and 
solve some puzzles of Nd systems |120j . Diversities are rewarding, but it is imperative 
to understand the relationship among various high precision potentials, and multitude of 
NN potentials has to be only a temporary situation. Will EFT give an answer? 
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The first challenge is to obtain full nuclear interactions at N3LO (derivation of the 
3NF in N3LO is in progress) and possibly even higher n. This implies understanding 
isospin violations and obtaining proper CD and CSB nuclear forces. Parameters in the 
EFT have to be uniquely determined and their physical meaning understood. EFT was 
successfully used in several FB studies |210|211j . some of relevance for BBN |212j . At- 
tempts were made to use EFT to improve Shell model [213]. More studies within EFT 
%PT are required. It is necessary to perform rigorous 3B, 4B... studies using NnLO and 
EM interactions. When we will have NnLO for n=4 and r3Bc, r4Bc etc done, then we 
will be able to conclude what are physical consequences when the calculations do not 
agree with the data. Therefore, a variety of several different approaches to calculate FB 
observables: momentum, configuration space, Faddeev and GFMC and NCSM should be 
encouraged. Relativistic effects were studied |214j at En = 28, 65, 135 and 250 MeV 
including relativistic kinematics, boost effects and Wigner spin rotation and as an input 
a relativistic NN on shell equivalent to AV18. The subject was reviewed at this confer- 
ence by F. Gross |215j . Relativity has to be fully and reliably taken into account at all 
levels: particle-particle and FB systems. The challenge is to obtain high precision baryon- 
baryon interaction. Progress is made by the Jiilich |216j and Nijmegen |217j groups and 
discussed at this conference |218|219j . FB studies are the ideal ground for searching for 
new physics. Examples are: i) studies of weak charge of the proton QW(p) = 1 — sin 2 Q w 
at JLab planned to achieve 0.3% accuracy of sin 2 w, (220j and ii) test of chiral symme- 
try breaking in systems with strangeness. The DEAR (DANE Exotic Atom Research) 
and SIDDHARTA (Silicon Drift Detector for Hadronic Atom Research by Timing Ap- 
plication) use precision X-ray spectroscopy of kaonic hydrogen and deuterium atoms to 
measure isospin dependent antikaon-N scattering lengths to test %SB |221j . FB research 
started with nuclear physics studies and - as a result of progress in theoretical studies and 
in technological development - a renaissance of nuclear physics is occurring. A special 
position of nuclear physics research at the intersection of particle, condensed matter and 
atomic physics and chemistry, biosciences and astrophysics makes it relevant now, when 
search for new physics is intensifying, when chemistry, atomic and nuclear physics employ 
ab initio procedures, when rare isotopes at the drip-line can be viewed as femtoscience 
and when astrophysics research critically depends on precise nuclear input. 1 So neutron 
superfluidity is relevant for phenomena occurring in the inner crust of the neutron stars. 
Precise 3 P 2 — Z F 2 pairing gap is important for cooling of neutron stars and cooling of young 
stars (T<10 5 y) is governed by neutrino emission (URCA process named after casino in 
Rio) |222|211j . Studies of neutron and proton drip lines, of Efimov effects and halo nuclei 
are a high priority and RIB will find applications in life sciences. 

Many topics, e.g. NN7 and pentaquarks, and a wide spectrum of atomic and chem- 
ical studies are left out. I just gave a partial list of challenges facing our FB research 
community. All tasks and challenges facing FB research merge into a comprehensive task 
of building the knowledge-based society, where science and physics in particular play a 
crucial role. Physicists should be like an English admiral in the anecdote describing the 
visit of Peter the Great to England. Peter was so impressed by the English Navy that 
he told the English admiral "If I were not the czar of all Russia, I would like to be an 
English admiral." to which the admiral replied "Sire, if I were not an English admiral, 
my greatest wish would be to become an English admiral." 
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